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Abstract: The color flux tube connecting a static quark-anti-quark pair in Yang-Mills 
theory supports massless transverse fluctuations, which are the Goldstone bosons of sponta- 
neously broken translation invariance. Just as in chiral perturbation theory, the dynamics 
of these Goldstone bosons is described by a systematic low-energy effective field theory. We 
use the effective theory to calculate the width of the fluctuating string at the 2-loop level, 
using both cylindrical and toroidal boundary conditions. At zero temperature, the string 
width diverges logarithmically with the quark-anti-quark distance r. On the other hand, 
at low but non-zero temperature T = 1/ j3, for r ^ f3 the string width diverges linearly. 
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1. Introduction 

Non-perturbative phenomena emerging in strongly coupled systems are quite common in 
physics. For example, in condensed matter physics, strongly coupled electrons are responsi- 
ble for quantum antiferromagnetism and high-temperature superconductivity, whose origin 
remains mysterious even after decades of intensive research. In particle physics, the con- 
finement of quarks is one of the most fundamental features of the strong interactions, and 
our understanding of the non-perturbative dynamics underlying quark confinement is still 
far from being complete. In general, it is very challenging to understand the dynamical 
mechanisms that are responsible for non-perturbative effects, also because collective modes 
play a crucial role. 

However, even if the origin of some non-perturbative effect remains unclear, one may 
be able to quantitatively understand some of its consequences. In particular, one may 
analyze the symmetries — internal or related to space-time — of the ground state of a 
non-perturbative system, and then work out the consequences of spontaneous symmetry 
breaking, even if the dynamics leading to this phenomenon may remain unclear. In this 
framework, the typical energy scale of the non-perturbative phenomenon represents an up- 
per bound for the energy of the phenomena that one wants to investigate. When a continu- 
ous symmetry breaks spontaneously, the resulting Goldstone bosons are the natural degrees 
of freedom to be taken into account in the low-energy regime of the theory. The Lagrangian 
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of the corresponding systematic low-energy effective theory consists of aU terms respecting 
the internal and space-time symmetries of the system. The terms in the Lagrangian are 
multiplied by low-energy parameters, whose explicit values can be derived only from the 
underlying quantum theory. They represent the high-energy non-perturbative input for 
the low-energy dynamics. 

The effective field theory approach is a very powerful tool for investigating strongly cou- 
pled systems. In particular, it provides a way to perform analytic calculations in physical 
systems that would otherwise be unmanageable. For example, chiral Lagrangians describe 
pion and pion-nucleon systems in a very accurate manner. Similarly, a systematic low- 
energy theory for magnons provides an excellent description of the low-energy dynamics of 
quantum antifcrromagnets. 

In 1980, Liischer, Symanzik, and Wcisz [1] proposed a low-energy effective string de- 
scription for the color flux tube connecting a static quark-anti-quark pair in the confined 
phase of d-dimensional Yang-Mills theory. During its time-evolution, the flux string sweeps 
out a 2-dimensional world-sheet, thereby spontaneously breaking the translation invariance 
in the transverse directions. The (d — 2) Goldstone bosons resulting from that breaking 
are the relevant degrees of freedom in the low-energy effective description of the string 
dynamics. The leading term of the effective action describes the dynamics of a thin string 
in the free-field approximation. Two non-trivial results follow from this observation. First, 
the leading correction to the linear term in the static potential [2] is — 7r(ci — 2)/(24r), where 
r is the distance between the two static sources. Second, the flux tube width is not fixed 
but increases like {d — 2) log(r/ro)/(27r(T) [3], with the string tension a and a length scale 
ro entering the effective theory as low-energy parameters. 

Many numerical lattice simulations have successfully demonstrated the validity of the 
predictions of the effective theory [4-19]. However, the very high numerical accuracy of the 
Monte Carlo simulations calls for improving the analytic results of the low-energy effective 
description. The sub-leading corrections to the free-field action of the {d — 2) Goldstone 
bosons depend on a number of parameters that, in general, can be fixed only by matching 
to the underlying Yang-Mills theory. Remarkably, Liischer and Weisz have pointed out that 
the low-energy effective string theory has a symmetry that had been overlooked before [20]. 
This new symmetry implies some relations between the values of the low-energy parameters 
of the effective theory. Interestingly, the Nambu-Goto action satisfies the constraints that 
one finds at next-to-leading order. In the special case of d = 3, the values of the low- 
energy parameters are completely fixed up to the next-to-ncxt-to-leading order and they 
turn out to be the same as those of the Nambu-Goto action [21]. In this paper we report 
on the computation of the correction to the increase of the flux tube width resulting from 
the sub-leading correction to the free-field action in d dimensions. The final expressions 
involve the Dedekind 77 function and the Eisenstein series E2 and E4. Due to the modular 
inversion property of those functions, we have the next-to-leading correction to the string 
width both at zero and at finite (but low) temperature. 

The paper is organized as follows. Section 2 contains the description of the systematic 
low-energy effective string theory for the dynamics of the confining string. In section 3 
toroidal boundary conditions, which describe a closed string wrapping around a compact 
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spatial dimension, are investigated. Similarly, section 4 discusses cylindrical boundary 
conditions, which correspond to the propagation of an open string that ends in the static 
quark- anti-quark charges. For both boundary conditions the width of the string is analyt- 
ically calculated at the 2-loop level. We present our conclusions in section 5. Technical 
details are described in four appendices. 



2. Effective string theory 

In this section we set up the framework for our computation. The Goldstone bosons 
resulting from the breaking of translation invariance are represented by a (d— 2)-component 
real-valued scalar field h(x, t) living in a 2-dimensional rectangular base-space (x, t) of size 
r X p. The field h describes the displacement of the vibrating string from the minimal 
length arrangement in the transverse {d — 2) dimensions. The leading-order free-string 
approximation is given by the following action containing two derivatives 

S2[h] = ^f^dt[ dxd^h-df,h, fi£{x,t}, (2.1) 
^ Jo Jo 

where a is the string tension. In the effective theory the string tension is a low-energy 
parameter whose explicit value can be obtained only from the underlying Yang-Mills theory. 
In our computation we consider toroidal boundary conditions, i.e. 

h{x,t) = h{x,t + f3); h{x,t) = h{x + r,t), (2.2) 

as well as cylindrical boundary conditions, i.e. 

h{x,t) = h{x,t + P); h{0,t) = h{r,t) = 6, (2.3) 

for the string. When one interprets the t-direction as Euclidean time, cylindrical boundary 
conditions describe the propagation of an open string at finite temperature. On the other 
hand, when one interprets the x-direction as Euclidean time, the same boundary conditions 
correspond to the propagation of a closed string that is created at "time" x = and 
annihilated at x = r. This dual interpretation gives rise to an open-closed string duality 
of the string theory. 

The first bulk correction to the free-string action contains four derivatives and is given 

by 



a [ dt [ dx \c2{d^,h ■ d^hf + C3{df,h ■ d^hf] . (2.4) 
Jo Jo 



The open-closed string duality [20] constrains the value of the two low-energy parameters 
C2 and C3 by 

, , ^ d - 4 , , 

(d - 2) C2 + C3 = -— . (2.5) 
o 

Note that, for d = 3, i.e. for a 1-component field h{x,t), the two terms in eq. ( p.4| ) are 
identical and hence the first correction is completely fixed. As shown in [21], for general 
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d a generalization of Liischer and Weisz's argument provides a further constraint on the 
next-to-leading order coefficients 

C2+C3 = -l. (2.6) 



Thus, for any value of d, the two independent constraints eq. (^) and eq. (|2.6|) completely 
fix the effective action at this perturbative order with C2 = | and C3 = Interestingly, 
the expansion of the Nambu-Goto action 

SNG[h] = a dt [ dxJl + dJi- dji + dth ■ djl + {d^h x dth)^ (2.7) 
Jo Jo 

satisfies these two constraints. 

Since cylindrical boundary conditions explicitly break translation invariance in the x- 
direction, one would expect surface terms (located at the boundaries at x = and x = r) 
to appear in the effective action. Remarkably, as was shown by Liischer and Weisz, due 
to open-closed string duality such terms are absent at leading order (i.e. there are no 
boundary terms with two derivatives). Boundary terms with an odd number of derivatives 
are excluded by parity symmetry. However, boundary terms with four derivatives do indeed 
exist. Fortunately, such terms contribute at one order higher than the four-derivative terms 
in the bulk that we discussed before. As a result, in our study boundary terms need not 
be taken into account. 

The squared width of the string is defined as the second moment of the field h, i.e. 

= Hh(..t) - T^f) ^ '^'"■!> - ,2,8) 
' " ^ ' ' ' fVh ej(p(-S[ft]) ^ ' 

Here S[h] is the effective string action and 

ho = ^f dt f dxh(x,t) (2.9) 
(3r Jo Jo 

is the equilibrium position of the string. For cylindrical boundary conditions, we have 
ho = 0- At next-to-leading order, the string action is given by S[h] = S2[h] + 5*4 [/i]. 
Similarly, at next-to- leading order, the field is replaced hy h{x, t) — > h{x, t)+ad^d^h{x, t), 
where a is a low-energy parameter. Expanding around the free-string action, the squared 
width of the string is given by 

w\x,t) = wl{x,t) - {{h{x,tf -hl)Si)o + 2a{{d^h{x,t) f)o 
+ {(.df,df,h{x,t) f)o - — j dtdx {ho ■ d^d^h{x,t))o 



j dt dx dt' dx' {df,d^h{x, t) ■ d^>d^>h{x', t') )o. (2.10) 



Here (• • • )o represents the vacuum expectation value with respect to the free-string action 
and 

wl{x) = {h{x,tf)o-{hl)o (2.11) 
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is the result for the squared width at leading order. 

Let us define G{x,t; x' ,t') = {h"'{x,t)h"'{x' ,t'))o as the free field propagator of a single 
component a. One then obtains 

{h{x, tf 54)0 = A{d-2){[{d- 2)c2 + C3] Ti + [2c2 + (d - 1)C3] T2 } . (2.12) 

Using the two constraints eq. (^^) and eq. ( |2.6D this implies 

{h{x,tf Si)o = ^^^Y^iTi - 2T2) -{d- 2)Ti. (2.13) 

The two terms Ti and T2 are given by 
r-r 

Ti = lim / dt' dx d^'G{x,t;x',t') df,>>G{x",t",x,t) d^>d^»G{x',t';x",t"), (2.14) 
e.e'^OJo Jo 

as well as 

T2= lim / dt' dx' df,'G{x,t;x',t') du"G{x",t";x,t) d^>du>>G{x',t',x",t"), (2.15) 
e.e'^OJo Jo 

where x" = x' + e and t" = t' + e' . Since they are ultraviolet divergent for {x',t') = {x",t"), 
the integrals defined above have been regularized using the point-splitting method. Finally, 
we have 

{h{x,t)-df,:df,>h{x',t'))o = {d-2)df,:df,:G{x,t;x',t'), 

{d^d^h{x, t) ■ d,>dA{x\ t'))o = {d-2) d,,d,,d,>d,>G{x, t- x', t'). (2.16) 



3. Toroidal boundary conditions 

In this section we present the computation of the string width with toroidal boundary 
conditions at next-to-leading order in the low-energy effective theory. As we show in 
Appendix the single-component free field propagator can be written as 

0(,.«) , ^ + -Lf-,« (j!I!H) + ^" ^ K, (3.1) 

^ ' ' 2aBr 2txg ^ \ r J n(l - g" ' ^ ^ 

n=l ^ ' \ J y 



« = g = e-2™, K = ^ —\og,^{iu). (3.2) 

r 12crr vrcr 



where 



It should be noted that t £ [0,(3]. In eq. ( |3.1| ) we have used translation invariance, i.e. 
{h{x,t) h{x' ,t')) = G{x — x',t — t'). Translation invariance also implies that the string 
width w{x) = w does not depend on the position x. At leading order, the squared width, 
wf^, is ultraviolet divergent and we regularize it using the point-splitting method 

wl= \ija{h{x,t)h{x',t')),-{^)^ = {d-2)[G{e,e')- f f G{x,t)\ . (3.3) 
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Using eq. ( B.2 ) one immediately obtains 



G'(e,e') = 7r^log-, T^ = 2'K^f^^'■ 



as well as 



such that 



P Jo r 



/5 



I2ar 



+ K, 



d- 2 , r 
-7, log — 



d-2 



log r]{iu). 



(3.4) 

(3.5) 
(3.6) 



The quantity rg is a low-energy parameter of dimension [length] . 

Let us now consider the corrections to this behavior resulting from the next-to-leading 
term, ^4, of the effective string action. By explicit calculation, it turns out that (/iq 54)0 = 
0. Hence, it remains to evaluate eq. ( 2.13| ) with Ti and T2 given by eq. ( p. 14 ) and eq. ( |2.15| ). 
We find 



1 



27rcr2 f3r 



log 2 log r](iu) 



1 



rj. ^Ti TTuE^jiu) _ E2{iu) 
^ 2 72a2r2 12cj2r2 ^ Ma^l3r' 



(3.7) 



Furthermore, using the two identities 



d^d^G{x -x',t- t') = -dtdtG{x -x',t- t') + 



dx'dx'dxdxG{x — x',t — t') = dt'dt'dtdtG{x — x',t — t') 



(3.8) 



we find that the terms proportional to a cancel and those proportional to vanish. Hence, 
at next-to leading order the squared width of the string is given by 



w 



1 



a (3r 



2 , ((i-2)2 /vr . ^ ,. ,,2 uE2{iu) , 1 



— [uE2iiu)Y 



+ 2^ 



(3.9) 



It should be noted that, order by order, this expression is modular invariant. In particular, 
it is invariant under the interchange of r and /3. 



4. Cylindrical boundary conditions 

Together with Appendix |d|, this section contains the calculation of the string width at 
next-to-leading order for cylindrical boundary conditions. In Appendix |^, we show that 
the single-component free field propagator can be written as 



with 



G{x,t;x' ,t') = — > sm sm ^— , (4.1) 

vrcT ^-^ \ r / \ T I nil — q^) 

n=l ^ ^ \ J. / 

u = ^, g = e-2-. (4.2) 
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In eq. (4J) we have used that {t — t') G [0,13]. We now calculate the string width w{r/2) 
at the midpoint x = r/2. It should be noted that, due to translation invariance in the 
i-direction, the string width does not depend on t. At leading order, the squared width 
wf^{r/2) is ultraviolet divergent and is again regularized using the point-splitting method. 
It turns out that 

wlir 2) = log - + log 4.3 

where the low-energy parameter ro is now given by 

vr 



ro = -Ve2 + e'2. (4.4) 

For j3 r, the second term on the right-hand side of eq. (|4.3| ) gives only exponentially small 
corrections to the leading logarithmic increase of the string width. The regime r ^ /3 [16] 
can be obtained using the inversion transformation rule given by eq. ( |A.13 ). Then we have 

„,?„W2) = |z^l„,A + il^. + 0(e-/.,. (4.5) 

Interestingly, this equation shows that at finite but low temperature, the squared string 
width increases linearly with the distance. Similar to the case of toroidal boundary condi- 
tions, we have to evaluate eq. ( p. 131 ). Since eq. 

(O) is well-defined only for {t -t') e [0, /3], 
it is convenient to split the integral over t, i.e. dt = dt + dt. We then obtain 

oo oo oo , 

^ ^ ^ ^ ^ ' \ ^ ^ f T^ke' /r I k TTke' /r 

ra=l m=l k=l 



^1 = ^ E E (-1)""" E Y^k + 'i'- 

n=l m=l k= ' 

r , ■Ki2kx + ke) r 

X dx cos 

Jo r L 



TT(2kx + ke) r 7r(2(m — n)x + (2m — l)e) 
cos 



-|- cos 



(l-g2n-l)(l_g2m-l) 

7r(2(m -\- n — l)x + {2m — l)e 



— r dt ^Q-2n{n+m-l)t/r ^ ^-2nin+m-l){(3-t)/r^ 



1 



(l-g2n-l)(l_g2m~l) 

as well as 



n=l m=l 
1 



(l-g2n-l)(l_g2m-l) 

E2{iu) 

gg^2j,2 



j (H ^g-27r(n-m,)t/r^2m~l _|_ g-27r(m-n)t/r^2n-l^ 



(4.7) 



The last term is the contribution of the integration dt after the limit e' — t- has been 
taken. The explicit evaluation of Ti and T2 is rather involved, and is thus relegated to 
Appendix |D[ 
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The two identities 



dxdr,G{x,t;x' ,t') = -dtdtG{x,t;x',t'), 
dx>d^idxdxG{x,t;x',t') = dt'dt'dtdtG{x,t;x' ,t' 



(4.8) 



imply that {{d^h{x,t))'^)Q = and {{d^d^h{x,t))'^)Q = 0. Hence the contribution to the 
squared width coming from the terms proportional to a and vanishes at next-to-leading 
order. Finally, inserting eq. ( p. 12 ) and eq. ( p.2| ) in eq. ( p. 13 ) we obtain 



w\t/2) = wl{r/2) + [E2{iu) - mi2iu)] (wl{r/2 



d-2 
Ana 



+ 



{d-2)-K 
12cj2r2 



d d-2 



12 



E2{iu) [E2{2iu) - E2{iu)] 



-Eoiiu 



,}. (4.9) 



5. Conclusions 



By now, a lot of numerical evidence supporting the validity of the low-energy effective string 
description of the long-distance static quark-anti-quark potential has been accumulated in 
lattice Yang-Mills theory. In several cases, the numerical data are so accurate that higher- 
order corrections to the leading free string approximation must be taken into account. 
Remarkably, open-closed string duality completely determines the terms in the effective 
action at next-to-leading order, without any additional low-energy parameters. In this 
paper we have presented the details of an analytic computation of the width of the color 
flux tube in the low-energy effective theory at next-to-leading order. Our result has been 
crucial for accurately describing the width of the color flux tube obtained in numerical 
simulations of Yang- Mills theory [18]. The results are expressed in closed form in terms of 
the Dedekind r/ function and of the Eisenstein series E2 and £^4. The modular inversion 
transformation property of those functions yields the next-to-leading order correction to 
the width both at zero and at finite (but low) temperatures. The calculation has been 
performed using both toroidal and cylindrical boundary conditions. 

The effective theory that we used describes string fluctuations in the continuum. In 
order to apply the results of our calculation to lattice field theories, one must be sufficiently 
close to the continuum limit. Before one reaches the continuum limit, the confining string 
in a lattice Yang-Mills theory is also affected by lattice artifacts. First of all, at very strong 
coupling the world-sheet swept out by the lattice string is rigid, i.e. it follows the discrete 
lattice steps and does not even have massless excitations. Only at weaker coupling, after 
crossing the roughening transition, the string world-sheet supports massless excitations 
and thus becomes rough. Consequently, the effective theory is applicable only in the rough 
phase. Since the lattice theory is invariant only under discrete rotations and not under the 
full Poincare group, before one reaches the continuum limit additional terms proportional 
to X]^=i 2(^M^M^)^ ^'^d X]/i=i 2('^/i^)^ enter the effective action in the bulk. Since these 
terms contain four derivatives, they are of sub-leading order. Hence, they have no effect 
on the Liischer term or on the leading logarithmic behavior of the string width. As a 
result, the Liischer term is completely universal. Provided its world-sheet is rough, even a 
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lattice string supports exactly massless modes which contribute — 7r/24r to the static quark 
potential. In order to incorporate lattice artifacts in the effective theory in a systematic 
manner, one must investigate whether additional boundary terms arise in the effective 
action. One must also reinvestigate the consequences of open-closed string duality, which 
were derived in [20] assuming full Poincare invariance. These are interesting problems for 
future studies, which may eventually be important for the correct description of numerical 
simulation data when the lattice spacing is not sufficiently small. 
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A. Infinite sums and products 

In this appendix we list the infinite sums and products which appear in our calculation. 
Some of them can be expressed in terms of the Dedekind rj function and the Eisenstein 
series. They are respectively defined by 

oo 

ry(r) = g 2^ n ' '? = e''"" (A.l) 

n=l 

and 



AU °° „2fc-l„n 
" n=l ^ 



(A.2) 



where are the Bernoulli numbers, defined through the expansion 

(2k)\ 



oo „ 



They are also related to the Riemann C, function of positive even integers and negative odd 
integers 

In the one-loop calculation one encounters the sum 

CXD _\ ^ OO OO 1^,y-^ OO 

s: fr^ - E s: V = - ^ ""^ (' - = - "(^^ ■ (^''^i 

n=l ^ ifc=ln=l fe=l 

where <^(r) is the Euler function, related to r] by 

r?(r) = gi/^V(r) • (A.6) 
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In two-loop calculations one finds sums of the type 

n=l ^ k=l^ ^ ' 

where the first identity is simply obtained by writing — q^) = X^fc^i'?"'' (hke in 

eq. ( |A.5| )) and then inverting the order of the two sums. The second identity follows from 
the definition of E2{t) or from eq. ( |A.5| ) using the relation 



TT = S^^W^ (A,8) 
?7(r) 12 



The identity (A. 7) also implies the following three equations: 

oo 2m— 1 1 

m=l ^ ' 



as well as 



and 



^(l + g^)2 - 8 + 6 24 ' ^^-^"^ 



The derivative of E2{t) can be expressed in terms of other Eisenstein series using 

lj-qE2{r) = ^ [E2{rf - E,{r)] . (A.12) 

The Dedekind rj function and the Eisenstein series E2k obey the following transformation 
rules under the inversion r — )• —1/t 

V{r) = ^v{—), E2{t) = \ E2{—) - — . (A.13) 

as well as 

^2fc(^) = E2k{^) , k>l . (A.14) 

These relations ensure fast convergence of the expressions in both the regimes /3 3> r 3> 
(the zero-temperature limit) and r ^ f3 ^ (the finite (but low) temperature 

case). 

B. Propagator on the torus 

The free field two-point function on a torus of size r x (3 satisfies 

-AG(x, t; x't') = -5{x - x') 6(t - t') ^ . (B.l) 

a a pr 
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Here the term \/a fir accounts for the zero- mode subtraction, which makes the Laplace 
operator A invertible in the orthogonal subspace. The solution of this equation can be 
expressed in closed form in various ways. We begin with the explicit expression given in 
Appendix A of [5] 

^ {t-f? 

- - ? — I - I 

r 



1 



2Txa 



fc=i 



with 



z = X - x' + i{t - i!) , Gk 



2k 



2a 



/3 



(B.2) 



E2k{iu) , u = - , 



(B.3) 



where the Eisenstein series E2k(j) are defined in our appendix^. Inserting this definition 
in eq. (B^) we get 



2-KaG{z) 



7r{Qmzy 



I3r 
2^e 



, 2-Kz ^Bk(2-Kf^ 1 



k=l 



{2k)l 2k \r 



2k 



2fc ^ ^2fc-l„n' 



n=l 



The terms enclosed in the first square brackets can be re-summed using eq. (A.3) 



xk^k' 

(21-,- -- 

k=l ^ ' n=\ 



(2fc)! 2k 



n 



w 

— +logw 



log(l 



— +logu;, 



(B.4) 



(B.5) 



with w = —2TTiz/r. Interchanging the order of the two sums in the last bracket of eq. (jBT 
we recognize the Taylor expansion of the cosine and arrive at the simpler expression 

TT{t-t'){t-t' - 



2TTaG{z) 



+ 3f?e 



(3r 
logf 1 — e 



oo —In 
2iTiz/r\ , , \ " Q 



^ ^ 1 - 



n=l 



27mz 

cos h 2 log (p{iu) 



(B.6) 



where the Euler function ip has already been defined in Appendix Using the identity 
( |A.5| ) and expanding the first logarithm in powers of e^'^*^/'', we obtain the Gaussian cor- 



relator on the torus in the compact form of eq. (3.1), where we have put G{x,t) = G{z). 
It should be noted that this expression converges in the range q <1 and < t — t' < /3. 



C. Propagator on the cylinder 

The Gaussian correlator G{x,t;x' ,t') = {h{x,t) h{x' ,t'))Q on a cylinder of size r x /3 with 
fixed boundary conditions at x = and x = r and periodic boundary conditions in t with 
period /3 can be conveniently written in terms of correlators G{z — z') and G{z + z'*) on a 
torus of size 2r x /3 (with z' = x' + it' and z'* = x' — it') 



oo 

1 . nnx 



G{x, t; x', t') = G{z - z') - G{z + z'*) 

oo 

E 

n=l 

This is the expression used in section 4 



imx' e-™(*-*')A 



+ e" 



-TTn{l3-t+t')/r 



sm sm 

vrcj ^ — ' r r 

n=l 



n(l 



(C.l) 
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D. Evaluation of Ti and T2 on the cylinder 



In this appendix we evaluate explicitly the terms Ti and T2 given in eq. ( ^^ ) and eq. • 
Performing the corresponding integrations one obtains 



^ 1 TT ^ ^ 

T2 - -Tl — -^rp, — > 

Z O T T ^ — ' 

m=l 



12 ^ ^(l-g2m-l)2 



1 ^2(^^) 

cT^r^ 96 



Using eq. ( |A.9| ), the sum may be written in closed form, which then yields 

1 ^ TT uE2{iu)[E2{2iu)-E2{iu)] 1 E2{iu) 

2 2 ^ cj2r2 144 cj2r2 96 

Performing the integrations in Ti and putting s = e"'^'^'/'", we obtain 



(D.l) 



(D.2) 



1 



1 22 



=1 



k I ^fc 



A; 7r(A; + 2n- l)e + 

cos — 



fc=ln=l 
,2/c+2n-l , ^2n-l 



+ 2n - 1 

„2fc+2n-l 



^2n-l 



1 



+ 



+ 



]^ _ (j2fc+2n— 1 ]^ _ (j2n 



A; + 2n-l vrfee s^+an-i _^ „fc+2n-i 2fc+2n-i 



— COS — 

k r 



^2n-l 



+ 



\ — qk+2n-l V 1 — g2fe+2n-l J. — q2n-l 
TT /3 ^ r(2m - l)g2"^-l(l +g2m-l) 



u2j,2 ^ 



m=l 



(1-9^ 



2m-l\3 



(D.3) 



The last sum may be written in terms of Eisenstein series using eq. ( [A. 11 ) or eq. ( |A.12 ). 
Introducing the two functions 



k 



k=l n=l 



k + 2n-l 



-a%' 



as well as 



Bia, 6) = f; f;(-l)^' + 2n - 1 ^„^,+2n^. 



k=l n=l 



and putting t = Q^'^i^+'>-^')/'^ ^ we can rewrite eq. ( p. 3D in the form 
1 ^(Ail,t) + A{s^t) 



1 99 



(D.4) 



(D.5) 



+ J2 [Mq^^Q^) + AiQ~^,Q^) + 2 A{q^,t) - B{q^,s q^) + B{q-^,s q^) + 2 A{1, q^ 



00 00 

+ 2Y,[Mq'^\Q')-B{q^,q^''') + B{q-^,q^+')]+2 J] A{q'^~^ ,q^)]] 



+ 



k=l 

TTU d 



k=l,k^j 



12cr2r2 dq 



q— [^2 (2m) - E2{iu)] . 



(D.6) 
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The two functions A and B can be written in terms of elementary functions 
aH 



A{a,b) 
B{a,b) 



(l-a2)(l + a6) (l-a2)2 
ab2 6(1 + b2) 

'(l-62)(l + a6) (1-62)2 



^^^ + ^')log(l + a6) + ^Mi:;^ + "^°^(^ + ') 



2(l + a)2 2(1 -a)2 



log(l + a6) . 



(D.7) 



The divergences of A{a,b) for a — )• 1 or 6 — )• 1 are only apparent. In particular, the useful 
ultraviolet limits to apply in eq. (p.6p are 



A{l,q) 



4(1 + . 



1 , 1 

T2 + «log- 



1 + . 



^(9,i) = 9- 



-g(l - (7) + (1 + log 



{l-q 



2\2 



^(I,t)+A(s2,t) 



2(1 
2r 



log' 



2r 



log 

16 8 Trie + ie 



7r|e + le' 
- + 0{e + ie' 



- + 0{e + ie') 



(D.8) 



We now rewrite this quantity in terms of Eisenstein series. We begin by explicitly writing 
the logarithmic terms of the double sum in eq. ( p.6| ), namely those arising from 



oo oo oo 

J]{2^[^((z^■+^<?^■)-i?(<?^<^^■+'=)+i?(g-^g^■+'=)] +2 Mq'-'^q')} , 



j=i k=i 

which are given by 



oo oo 



j=l k=l ^ ^ ' 



.j+k ^j+k 
log(l - q'') + 



+ 



+ 



(l + g.-^)2l°S(l-'^ 

log(l + 

k=l,kf^j 
oo 



(1 - q^+k)^ 



log(l + /) 



(1 _ g2(i+fc))2 

g'=(l + g2(fe-i)) 
qk^l _ qj-ky - ^ qj (I - q2(k-3))2 



log(l + q"") 
log(l + Q'=)]} 



J^logy 

k=l 

oo 



"~ ^ ^ _ qK 

\ (1 + q^Y ^ iTg^ ^ (1 + q^Y ^"^ 1 + q^ 



+ E 

CO 

E 



log- 



(D.IO) 



The last sum cancels exactly against the logarithmic terms resulting from the single sums 
(i.e. the second line) in eq. ( p.6| ). Apart from the first term in eq. ( P.IOD , the only remammg 
terms are those associated with eq. (p.8|). Putting all these terms together we obtain 



Ti 



V2 I (l + gi)2 ^ 8 



log 



2r 



Trie + ie 



(1+^ 

k\2 



(D.ll) 
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By comparison with eq. p. 3D in the second factor of the first term we recognize the contri- 
bution to the squared width of the flux tube at leading order. Applying the identity ( A.1C| ) 



and the definition (A.l) of the Dedekind rj function, we finally obtain 



if / rrP'iliu) l\ d 

^2[m^u)-^E,i2^u)] ^log - - j + 2vrng- [ii;,(2.n) - ^.M] 

(D.12) 



where ro has been defined in eq. ([4.4|). Inserting eq. ( D.12 ) in eq. ( [D.2 ) one immediately 



obtains an explicit expression for T2 in terms of Eisenstein series. 
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